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Abstract 


In this work we perform discrete modeling and analysis of the deformation of glassy 
amorphous polymer. Amorphous polymers consist of long molecular chains that are 
randomly oriented. Modeling of these chains is done by using a coarse grained self 
avoiding random walk (SAW) model. A chain generated by SAW can not cross itself but 
two chains can pass through a common position in space called the entanglement point. 
This simulates the effect of entanglements on the deformation. Care is taken to avoid a 
particular chain from overlapping another chain. Chains are obtained in a regular grid of 
dimensions obtained from specifications of the polymer under study. 

The final polymer model is like a nonlinear truss and the analysis of this structure takes 
account of material and geometric non-linearities. The constitutive response of the 
members of the truss is like rubber elasticity obeying Langevin statistics. Total 
Lagrangian formulation is applied to the model which uses the Second Piola-Kirchoff 
stress and the Green Lagrange strain tensor. 

Finite element analysis of the truss structure is done to obtain stress-strain curves. 

This model is capable of providing quantitative information on 

(1) stress-strain behavior, molecular extension , chain fracture and slippage 

(2) importance of the molecular weight and its distribution 

(3) effects of testing variable, such as temperature and strain rate. 
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CHAPTER 1; INTRODUCTION 


Now a days polymers are widely used in various fields, such as structures, composites 
etc. The most important aspect of polymer material is the diversity of accessible 
molecular composition, as well as a choice of monomers. Establishing a good connection 
between molecular composition and material properties is a difficult task due to complex 
structural and dynamic features observed in these materials. Our work will be based on 
predicting a model and simulating it to get the properties of polymer under study. In 
particular it is concentrated to get stress-strain curve to predict the mechanical properties 
of amorphous polymers. 

Polymers consist of large molecules i.e. macromolecules. A monomer is the repeating 
unit of which a polymer is made of, for e.g. propylene is the monomer which polymerizes 
to give polypropylene. Polymer crystals show direction-dependent (anisotropic) 
properties. For e.g. Young’s modulus of polyethylene in chain direction at room 
temperature is 100 times of the Young’s modulus in transverse direction, this difference 
in modulus is due to the presence of two types of bond present in the crystal: 1) strong 
and stiff bond along chain axis and 2) weak and soft secondary bond acting in the 
transverse direction. However these weak (Van der Waals) bonds are also present with 
strong bonds along chain axes. 

When a polymer is cooled down from molten state it either crystallizes or cools down to 
its glassy amorphous state as shown in Figure 1. The temperature at which the slope in 
volume-temperature graphs changes is referred as glass transition temperature, . 

Basically the chains in amorphous polymers are randomly oriented and that of crystalline 
polymers are oriented in a regular fashion. In our work we had used the kinetic model 
developed by Termonia and Smith. [1,2] for tensile deformation of amorphous polymers. 
This model is based on Eyring Activation rate theory [3] and classical theory of rubber 
elasticity [4]. 
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Figure 1 Cooling of molten polymer 


1,1 Comparison with other models: 

The macroscopic or bulk models take into account the overdl stress-strain behavior of a 
polymer and not local stress-strain variation which occurs at microscopic or atomistic 
levels. On the other hand atomistic models consider individual atoms and bonds 
explicitly which makes the study of large system almost impossible. The model which we 
are using considers the local extension ratios and stresses of chain strands by applying the 
concepts of rubber elasticity theory to each individual rubbery chain strands and linear 
elasticity theory to weak Van der Waals bonds. The chain segment between two 
entanglement points contain hundreds of atoms in rubbery phase, therefore the model is 
neither macro nor atomistic but is mesoscopic. This model a) includes the most relevant 
details in simulations of mechanical properties of polymer and b) keeps the computer 
time and memory requirements within acceptable limits. 

The model is capable of giving information on 

1) Stress-Strain behavior based on molecular extension, chain fracture and slippage. 

2) Effect of testing variables, such as temperature and strain rate. 
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3) Importance of molecular weight and its distribution. 


1.2 Literature Survey: 

Our model is based on the work done by Termonia and Smith [1]. Many models have 
been developed to describe phenomenon like brittle fracture, yielding, strain hardening 
etc [4]. This model is based on Eyring activation rate theory [3] which takes into account 
weak attractive forces between chains and also chain slippage through entanglements. 
According to our model the weak Van der Waals bonds break according to the kinetic 
theory of fracture [5]. The stresses in the strands are obtained using the classical theory of 
rubber elasticity [4]. The simulation of Van der Waals bond breaking is performed with 
the help of Monte Carlo lottery procedure developed by Termonia [1]. 

The theory of Rubber Elasticity was developed by Meyer, Von Susich and Valko [4] and 
Karrer [4]. Karrer developed the concept of Brownian motion at various portion of long 
chain molecule structure. Meyer, Von Susich and Valko considered the relation between 
elasticity and the ability of chain to take up an irregular or statistically determined form 
by virtue of the interchange of energy between its constitutive elements and surrounding 
atom. This is fundamental concept of now theory of rubber elasticity. 

Busse [6] considered the general condition necessary for rubber like elasticity and gave 
the concept that rubber like elasticity is associated with micro-Brownian motion of 
individual element of long chain molecules. Later on this concept helped to develop the 
concept of Van der Waals bonds in polymer structure. 

When the kinetic theory of elasticity was first proposed Meyer and Ferri [4] realized the 
importance of variation of elastic tension with temperature i.e. deformation is entropy 
driven. In the absence of external forces the chain molecules will adopt configurations 
corresponding to a state of maximum entropy (contracted configuration of chain). When 
forces are applied the chains will tend to extend in the direction of the force, so reducing 
the entropy and produce a state of strain. Our model is based on random walk, and the 
solution to random walk (chain) problem was first given by Kuhn [7] and by Guth and 
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Mark [8]. The problem consists of essentially in evaluating the relative number of 
configuration of the chain, for given values of x, y and z or getting the information on the 
probability of a chain having a particular end-to-end length. 

Statistical theory of a polymer network is first proposed by Kuhn [7] who showed that the 
elastic modulus could be directly related to the number of molecules per unit volume of 
material. Later on Wall [4], James and Guth [13], James [10], Treloar [11] made their 
contribution to statistical theory of rubber elasticity which was based on Gaussian 
statistics. Later on non-Gaussian Statistical network theory was developed which takes 
into account the finite extensibility of the chain and hence of the network. Kuhn and Grun 
[12] gave the idea of statistical treatment of random chain by non-Gaussian statistical 
network theory. They used the inverse Langevin function to get an expression for 
probability of a chain having a particular end-to-end distance. James and Guth [13] used 
the inverse Langevin distribution function to give expression for the force per unit area. 

1.3 Overview of the Work: 

Present work is an attempt at predicting the mechanical behavior of glassy amorphous 
polymers from a mesoscopic perspective. Modeling of amorphous polymer chains is done 
by using a coarse grained self avoiding random walk (SAW) model. The final polymer 
model is like a nonlinear truss and the analysis of this structure takes account of material 
and geometric non-linearities. Total Lagrangian formulation of truss is done based on the 
Second Piola-Kirchoff stress and the Green Lagrange strain tensor. The polymer model is 
assumed to behave according to the laws of rubber elasticity. Finally strain energy-strain 
and stress-strain curves were obtained for the polymer model. 
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CHAPTER 2: RUBBER ELASTICITY 


2.1 Introduction: 

In this work we intend to build a lattice based, mesoscopic model of deformation. This 
model will be based on the statistical theory of rubber elasticity. Analogy will be drawn 
between the affine deformation of networks of rubber macromolecules and plastic 
deformation of amorphous polymers. This analogy is motivated by the work of Arruda 
and Boyce [14] and Wu and Van der Giessen [15] on constitutive models of glassy 
polymers. 

2.2 Statistical properties of long chain molecules: 

The Theory of Rubber Elasticity on the basis of structural and thermodynamic concepts is 
developed in two stages: 

1) Development of mathematical expressions for probability and entropy as function 
of end to end distance. 

2) Treating network of such long chain molecules by statistical methods. 

Ensemble of random chains can only be described by means of a spatial distribution 
function. A random chain is shown in Figure 1. 



Figure 2 Random chain 
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The end-to-end distance(r) and the radius of gyration(s) are two different measures of 
random chain. The former is the distance between the chain ends and the latter is the root 
mean square distance of the collection of atoms from there common center of gravity. 
The following relation holds good between r and 5 . 

<s^>=<r^>l6 ( 1 ) 

2.3 End-to-end distance of a macromolecular chain: 

Average end-to-end distances of long macromolecular chains can be calculated using 
various levels of approximations. A polymer molecule in a ‘theta’ solvent (solvent that 
does not cause the coil to either expand or shrink) follows the equation: 

<r^>o=Cnl^ ( 2 ) 

Where the subscript ‘0’ signifies the theta condition, n is the number of chain links in a 
particular chain and 1 is the length of each link and < > signifies average over the value. 
The constant G depends on the nature of the polymer. This formula considers only the 
short range interactions. If the long range interaction are considered then, 

<r^ >=a^ <r^ >Q, (3) 

where a is the expansion factor and is affected by temperature and type of solvent. 

2.4 Models for calculation of average end-to-end distances: 

Some models for calculation of end-to-end distance of molecular chain are as follows 

1) Freely jointed chain: 

For a freely jointed chain with main chain bonds (segments) each bond having length 1, it 
can be shown that, 

< >= nf . (4) 

In this approximation it is assumed that the orientation of the individual segments is 
totally uncorrelated. See Figure 2. 
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Figure 3 Freely jointed chain 


2) Freely rotating chain: 

A freely rotating chain model assumes that the bond angle ^ is a constant. No particular 
chain conformation is preferred and the average projection of bond i+1 along a direction 
perpendicular to bond i is zero. Taking 6 « 1 10" this assumption yields, 

<r^>^2nl\ (5) 

It is assumed that each C-C bond permits of free rotation, so that, for example, the bond 
CjCj may be considered to rotate about as axis, while C 3 C 4 rotate about CjC^ as 

axis, and so on. 



Figure 4 Freely rotating chain 
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3) Hindered rotation chains; 

Above models neglect any hindrances to rotation. They assume that all positions of 
rotating bond are equally possible, it doesn’t respond to conditions in actual molecular 
chain where attractive and repulsive forces between atoms of neighboring chains makes 
certain positions of rotation more probable. These effects are conveniently described in 
terms of difference in potential energy between different positions. 


For example ethane ( CH-^ CH ^ ), where hydrogen atoms repel each other causing energy 
maxima in eclipsed position and minima in staggered position. For this model 

( 6 ) 


<r 


2 ^^^ 2 r l + cos(;r-<9) l+<cos^ > 

1-cos(;t— 0) l-<cos(z)> 


For freely rotating chain model, <cos^ >= 0. For higher alkanes the quantity <cos^ > 
can be obtained from consideration of the existence of trans and gauche states. <cos^ > 

is function of potential energy and temperature , as temperature increases potential 
becomes less effective and RMS length tends to approach that of unrestricted rotation. 



Staggered position 
(most stable) 




(least stable) 


Intermediate position, 
definition of torsion angle(^ ) 


Figure 4 Rotational isomers of ethane from a view along the C-C bond: carbon-shaded; hydrogen- 
white 
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4) Concept of an equivalent chain: 

The constant C in equation (2) arises from short range interactions. We can actually 
assume an equivalent chain comprising of N hypothetical segments each of length/ , so 
that the equivalent chain behaves like a freely —jointed chain with 

< >= N'P (7) 

and<r„^>=W/'. (8) 

Thus for real polymer with r^^=0.83Nl and<r' >o=6.7M' , with N/N'^\0, this 

means that each segment of the equivalent chain consists of 10 real bonds. In other 
words, the actual chain orientation becomes uncorrelated after every 10 bonds. 



Figure 5 Equivalent chain 


2.5 Distribution of end-to-end vectors: 


The distribution of end-to-end vectors for Gaussian chains follows the distribution 


function: 


P(x, y, z)dxdydz 


371 

zn 

3(x^ +/ +z^) 

>0 ^ 

CAp 

2<P>o 


dxdydz 


(9) 


Where < P >q is the average end-to-end distance of the chain. The above formula gives 

the probability of finding the chain, whose origin is at (0, 0, 0) at the point (x, y, z). This 
formula is strictly true for any chain that can be represented as a hypothetical equivalent 
chain with n freely jointed chain, each having length (1) as explained above. 
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2.6 Rubber Elastic State: 

2.6.1 Introduction: 

Unlike in metals elasticity exhibited by rubbers is entropy-driven. This is demonstrated 
by the thermo-elastic effects: a) stretched rubber subjected to uniaxial load contracts 
reversibly on heating; and b) rubber gives out heat reversibly when stretched. This 
confirms the fact that entropy of rubber decreases on stretching. This behavior is different 
from energy driven or energy elasticity of metals, ceramics. 

2.6.2 Thermodynamic consideration: 

To illustrate the contrast between energy and entropy elasticity thermodynamic relations 
are set for stretching of a sample uniaxialy. According to the first law of thermodynamics 
the energy change, dE, for a body subjected to a transfer of heat, dQ, and change in 
volume, dV, under pressure P and uniaxial deformation is given by 

dE = dQ-PdV+fdL, (10) 

where f is the retractive force and dl is the change in length, 1, of the sample. For the 

reversible case TdS = dQ, and so 

dE = TdS-PdV+fdL. (11) 

We define Gibbs free energy (G) as: 

G =H-TS, (12) 

Where the enthalpy H is given by, 

H=E+PV. (13) 

Hence 

dG=dE+PdV+ VdP-TdS-SdT. (14) 

From (10) and (13) 

dG=fdL+VdP-SdT. (15) 

Thus we have 

'—1 =/ and (16a) 

I dL 
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Ur A. 


=-5, 


(16b) 


Since the order of differentiation is unimportant: 


d 


dT 


\ A 


V 




dL 


ST), 


p Jpj 


implying 




ydLj 


P.T 


Partial derivative of G with respect to L at constant P and constant T is given by: 


\dLj 


P,T 


^dL 


-T 


Ax 


\SL XpT 


Combining equations (18) and (19) gives: 


/ = 


dL 


+ 


Ax 


vary 


pj. 


The first term in rhs of the above equation can be written as 


'm'] 

dL 


Ax 


V ST j pj 


\dLj 


PX 


Now assuming that volume is constant during deformation i.e. 

^ dV^ 


ydLj 


! 0 , we have 


px 


dH 

dL 


^ dE^ 


pX 


ydLj 


px 


And so we have equation for force as 

^ dE^ 


f = 


\SLjp.p 


+T 


ydTj 


pJ. 


(17) 


(18) 


(19) 


(20) 


( 21 ) 


( 22 ) 


(23) 


The first term in the equation (23) is associated with change in internal energy 
accompanying deformation at constant pressure while the second term accounts for the 
change in entropy by deformation as: 
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(24) 


from (18). 


UrJ 


p,L 


^dL 


P,T 


So equation (23) shows that the total force has an energetic and an entropic term. The 
entropic part, in case of polymers arises from the changes in chain orientations 
accompanying deformation. 

2.6.3 The Statistical Mechanical Theory of Rubber Elasticity: 

To develop the concept of Statistical Mechanical Theory of Rubber Elasticity we have to 
go through the concept of affine and phantom models. In affine deformation model the 
junction points are assumed to have fixed positions defined by deformation ratio, while 
the chains between junction points are free to take any possible conformation. The 
junction points of phantom network are allowed to fluctuate about their mean position 
and the chains between these junction points can take great many possible conformations. 
Statistical Mechanical theory starts from affine network model which obeys following 
assumptions; 

1) The chain segments between crosslinks are represented by Gaussian statistics of 
phantom chains. 

2) The entropy of the network is the sum of the entropies of the individual chains 
which are N in number per unit volume. 

3) The deformations on the molecular level are the same as that on a macroscopic 
level, i.e., deformation is affine. 

4) The. unstressed network is isotropic. 

5) The volume remains constant during deformation. 

The distribution of end-to-end vectors is as follows from (9) 


P{x, y, z)dxdydz 



312 

3(x^+y+z')’ 

^2;r<r^ >oj 

CAp 

2<r^>o 


dxdydz . 


Entropy is given by the equation 

S = k\nP , 


(25) 
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and so 


S = k^-ln 

2 1 2;r < >„ 


+z^) 
2<r' >„ 


Therefore 


S = C-k 


2<r^ >/ 


where C is a constant; k is the Boltzmann’s constant. 

An unstressed chain is represented by end-to-end vector starting at the origin and having 
position vector rQ={x^,yQ,Zf^)mA end-to-end vector /* = (;c, z) corresponds for 

position vector of stressed state and is related to through relations; 

X = / 1 |Xq y = A2T0 ^ ~ ^^0 

The change in entropy of the chain due to this stretch is given by 


AS^ ='^AS = -3k 




2<r’>. 


Here the change in entropy for the network is the sum of the contributions of all the chain 
segments of the network. 

It is assumed that original system is isotropic i.e. 


N N 




N N 


E^o+Z^o+Z^o =Z^o 


N N N 1 ^ 

I I 1 j 1 


N<r^ >„ 


where N is the number of Gaussian chain segments in the system. From equation (27) and 
(28) we have entropy change for network as; 

AS,=-^kT(A,^+Ji^+A^-3). 

Consequently the change in internal energy is given by 


AG = -TAS^ = -j-NkT(2^+jq -3) . 
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Assuming volume to be constant during deformation we get: 

W=l- 

For uniaxial stress deformation with 

>^ = 1 , we have from (29) 

^ = /^ = , where A=—. 

V2. Lq 

Here L is the current length in x, direction in which the original length was Lq . 
Using equation (16a) we get 


( 29 ) 


/ = 


V Jry 


dX dL 


implying 


/ = 


NkT 


A- 


1 


(30) 


Cauchy stress is obtained using A = Aq/ X (Aisthe area of cross section of specimen): 


<7 = ^ = nRT(x^-- 
A) V ^ 


(31) 


where n is the number of moles of Gaussian chain segment and Fg is the volume of 
specimen. Also n! Vg can be expressed as: 


n _ nM^ 1 


V, Vg M, M, 

where is the molecular weight of each Gaussian segment. Hence 


cr = 






(32) 


The rubber modulus is linearly dependent on temperature and inversely 


proportional to. 
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2.6.4 Drawbacks of classical statistical theory and implementation of non-Gaussian 
statistical theory: 

Gaussian statistics is good approximation only if end-to-end distance is considerably 
smaller than chain length (NI). Non-Gaussian statistical treatment takes into account the 
finite extensibility of the chain, and thus leads to a more realistic form of distribution 
function which is valid over the whole range of r values up to the fully extended length. 
Also Gaussian approximation becomes increasingly inadequate with increasing crosslink 
density. For very short chains, for example < 5 , the mean chain extension in the 
unstrained state already exceeds that for which the Gaussian approximation is valid. 
Therefore non-Gaussian treatment is necessary for accurate determination of behavior 
even at smallest strains. Such a treatment is adopted by Kuhn [7] and he showed that the 
probability density P(r) is given by 

In P(r) = const.- N —/3 + \n — - — , (33) 

_nl sinhy^J 



The complete force -extension relation is valid from r/Nl =0 to r/Nl =1, where the inverse 
Langevin function tends to infinity. In this work we will use the Non-Gaussian theory to 
characterize our network of polymer chains in the simulation. The Inverse Langevin 
function is approximated by the relation 

r'(x) = x(3-x')/(l-x'), (36) 

where x is any variable. Figure 5 shows the graph for Inverse Langevin and its 
approximation. The graph shows that as x — > 1 the function L ’ (x) — > . This behavior 
is shown by the strain energy-strain curve of polymer model. 
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Figure 6 Inverse Langevin function exact and approximation 
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CHAPTER 3: SELF AVOIDING RANDOM WALK (S.A.W) 


3.1 Introduction: 

A polymer molecule can take many different shapes (conformations). This is due to its 
freedom of rotation about C-C bonds; for example ethane has two conformational states 
called eclipsed (energy maxima) state and staggered (energy minima) state as shown in 
Figure 1. 



Staggered position Eclipsed position Intermediate position, 


(most stable) (least stable) definition of torsion angle( ^ ) 

Figure 7 Rotational isomers of ethane from a view along the C-C bond: carbon-shaded; hydrogen- 
white 


These conformation states are defined by the torsion angle between C-C bonds. Higher 
alkanes like n-butane have stable conformation states referred as trans(T) and gauche (G 
and G ). Higher alkanes exist in conformation states which are combination of these 
trans(T) and gauche(G and G ) states. Flory developed Rotational Isomeric State 
Approximation method, which is a convenient way for dealing with conformational states 
of polymers. In this method each molecule is treated as existing only in state 
corresponding to potential energy minima, i.e. different combination of trans(T) and 
gauche(G and G’) . Polymer chains exhibit in several cases a random (Gaussian) 
conformation, i.e. a random distribution of trans and gauche states. Random 
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macromolecular chains are found in solutions of polymers in good solvents, in polymer 
melts and also in glassy amorphous polymers. 

An ensemble of polymer chains can only be described by means of spatial distribution 
functions. Two different measures of random chains are commonly used a) end-to-end 
distance(rj and b) radius of gyration (s). Former is distance between chain ends and latter 
is root-mean-square distance of collection of atoms from their center of gravity. Hence 5 
is defined by 


is the vector from center of gravity to atom 1 of the polymer, m. is the mass of atom. 
For long chain polymers following Gaussian statistics the relation holds, 

<5^>=<r^>/6. (2) 


3.2 S.A.W: A sensible model for real polymers; 

Real polymer molecules exist in continuous space having tetrahedral (109.47°) bond 
angles, non trivial energy surface for bond rotational angles, and complicated monomer- 
monomer interaction potential, while S.A.W exist in discrete lattice with non-tetrahedral 
bond angles (e.g. 90°and 180° on simple cubic lattice), energy independent of bond 
rotational angles, and repulsive monomer-monomer potential. Inspite of these factors 
S.A.W is a perfect model for some aspects of behavior of linear polymers in good 
solvents (where the temperature of solvent is greater than theta solvents or solvents where 
coils expand). This behavior of S.A.W chain model arises from Universality, which plays 
a central role in modern theory of critical phenomenon. In short, critical statistical - 
mechanical systems are divided into small number of universality classes, which are 
typically characterized by spatial dimensionality, symmetries and other general 
properties. Near these critical points the leading asymptotic behavior is exactly same for 
all system of universality classes, details of chemical stmcture; interaction energies etc. 
are completely irrelevant. For this reason it is essential to determine universal and non- 
universal quantities. 
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To determine universal quantities, any mathematical model belonging to the same 
universality class as the system under study is taken, and solving it to determine the exact 
values of universal quantities. 

The behavior of polymer molecule as the chain length tends to infinity is a critical 
phenomenon in above sense. It is found that mean-square radius of gyration <s^ > of 
linear polymer molecule consists of monomers units has leading asymptotic behavior, 

<s^>=AN^''fiN). (3) 

Here f(N) is some function of N; as N -><x> where critical exponent v« 0.588 is 
universal, i.e. exactly the same for all polymers, solvents and temperature (temperature 
greater than theta temperature). The critical amplitude A is non-universal, i.e. it depends 
on the polymer solvent and temperature. Therefore linear polymer chains which exhibit 
some flexibility and has short range predominantly repulsive monomer-monomer 
interaction lies in the same universal class as that of S.A.W. 


3.3 S.A.W: notation and terms. 

Let L be some regular d-dimension lattice. A 7V-step (S.A.W) co on L is a sequence of 
distant point (0^,(0^, in Z such that each point is the nearest neighbor of its 

predecessor. N is the number of steps in (o . Following are the measuring terms of the size 
of N-step S.A.W. 

• The squared end-to-end distance of S.A.W. 

. (4) 


The squared radius of gyration 


5 =- 


1 


■E(6>,-^»y)'- 


2(iv+i)^tr 

The mean-square distance of a monomer from the endpoint 


(5) 


r? = — — yr(£u,-tuo)^ +(<y, • 

” 2(W+l)ts'^^ ' ^ -I 


( 6 ) 


The mean values <r^ >^,<s^ >^and<rl>^^ are believed to have asymptotic 
behavior 


19 





as TV-)- 00 ,where v is another (universal) critical exponent. Moreover the amplitude 
ratios; 


A 

_ 2 


<r > 


N 


(7a) 


R -HkllL 
~ 2 ^ 
<r 


are expected to approach universal values in the limit TV -» oo . 


(7b) 


3.4 Generation of the model of polymer solid: 

3.4.1 Introduction; 

The fundamental concept of rubber elasticity put forward by Meyer, von Susich and 
Voiko [4] developed the relation between elasticity and the ability of the chain to take up 
an irregular or statistically determined form by virtue of the interchange of energy 
between its constituent elements and surrounding atoms. Later on Karrer [4] noted the 
properties of long chain structures and found that in long chain molecules in which 
internal rotation about bonds can take place will be subjected to Brownian motion at 
various portions of its structure. This so-called micro-Brownian motion of the individual 
elements of long chain molecules is similar to the thermal motion of molecules of an 
ordinary liquid, with the single difference that the moving units, instead of being entirely 
independent, are connected together in form of chain. It is known that properties of the 
forces exerted by one molecule on another (secondary forces which are weak Van der 
Waals forces) are of the same order of magnitude as the intermolecular forces in an 
ordinary low-molecular weight liquid. These Van der Waals forces are necessary if the 
random rotations and the resulting changes in configuration of the molecules upon which 
their elasticity depends are to take place. 

3.4.2 Model Generation: 

The undeformed polymer solid is represented by a dense network of coiled, entangled 
macromolecules as shown in Figure 2. We start with a regular array of entanglement 
points depicted in Figure 3. Entanglement points axe topologically defined points in 
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polymer model through which a polymer molecular chain passes when it is in motion. 
The chain segments between entanglements are tied together by Van der Waals force or 
bonds (dotted lines). These secondary forces provide the necessary initial stiffness to the 
material. For convenience the three-dimensional network is given a planar conformation. 



Figure 8 Schematic representation of the undeformed polymer solid 



Figure 9 A more schematic representation of the network, in which the details of chain 
configuration are omitted altogether 
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Van der Waals bonds are represented by dotted lines and solid lines denote chain vectors 
between entanglements. In three dimensional space chain vectors are randomly oriented 
in an actual, undefoimed specimen i.e. < cos^ 9 >= 1 /3 , where 6 is the vector angle with 

the draw axis. In the two-dimensional representation we tzike the same angle as ^=54.7" 
for all vector orientation along the y-axis. 

3.4.3 S.A.W representation of a single chain: 

A single chain is represented by S.A.W; the algorithm used for generation of S.A.W is as 
follows (in our algorithm we start the chain from first entanglement point of the regular 
grid): 

Title: S.A.W single chain 
Subroutine: SC (N) 

Comment: This subroutine returns S.A.W chain 

start: for i=l to N do 

random nearest neighbor of 

if (o^ € {co} 

(D^=(D,_, 

goto start 
else 

end do. 

Here ( 2 >is the sequence of distant points o)^,co^, co^ where each member is the nearest 

neighbor of its predecessor, N is the number of steps of S.A.W (chain), co^ is he start 
point of the chain. S.A.W of 100 steps generated by using this algorithm is shown below 
in Figure 4. 
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Figure 2 Multiple chains generated 
3.4.4 Multiple chains generation: 

Multiple chains are generated by using the algorithm shown below: 

Title: Multiple Chain 

function: MC(fl) ) 

Comment: this routine returns multiple S.A.W chains 
start: call subroutine SC(N) 
if 
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present node and the node to which the chain is to move is in some 
previous chain then move to the nearest neighbor of present node which is 
not in any of the previous chains, 
end if. 

This algorithm takes care of the restriction on the chains that they should not 
overlap other chains. For modeling of polymer this restriction is required because 
no two chain strands can occupy the same position in space due to repulsion of 
the bonds. See Figure 5 

3.4.5 Applying boundary conditions to the polymer model: 

We had considered the polymer model to be taken out of a large polymer 
specimen which is continuous in x-direction and so we had applied periodic 
boundary condition in x-direction. However we have chopped the chains which 
are going out of the grid from extreme v-coordinates (See Figure 6). 



Figure 3 Multiple chains showing periodic boundary conditions applied to the model 

From the Figure it can be seen that the lines between nodes at extreme x- 
coordinates of the grid depicts the continuity of the grid patch in x-direction. 

3.4.6 Final model; 

Final model of the polymer is obtained by taking into consideration both Van der 
Waals bonds and the primary bonds which are due to C-C bonds. Here it is worth 
mentioning that the chain segments between two entanglement points (we call 
them strands) are made up of large number of atoms which constitute the polymer 
molecule. See Figure 6. In the figure the red lines show the Van der Waals bonds 
and the blue lines show the molecular chain obtained from S. A.W algorithm. 
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Figure 4 Final model of the polymer with Van der Waais bond and molecular bonds present 
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CHAPTER 4: SIMULATING SAW MODEL 


4. 1 Introduction: 

The model of polymer developed using S.A.W is like a truss structure, each element of 
this structure is like a spring that experiences large deformations and rotation, and the 
stress-strain relation of this element is also nonlinear, this give rise to geometric 
nonlinearity and material nonlinearity respectively. Hence when we analyzed the model 
using Finite Element Methods geometric nonlinearity and material nonlinearity were 
implemented in the analysis methods of the model. 

4.2 Finite Element Formulation: 

The S.A.W model is analyzed using corotational Lagrangian meshes [17] where the 
nodes and elements move with the material. Boundaries and interfaces treatment is 
simplified as they remain coincident with element edges. Lagrangian meshes are 
advantageous for history-dependent materials because in Lagrangian formulation 
constitutive equations are always evaluated at the same material points. 

Finite element discretisation with Lagrangian meshes are . classified as 

1) Updated Lagrangian Formulation and 

2) Total Lagrangian Formulation. 

Here both formulations use Lagrangian description, i.e. dependent variables are function 
of material (Lagrangian) coordinate and time. We are using total Lagrangian formulation, 
where the weak form involves integrals over the initial (reference) configuration and the 
derivatives are taken with respect to material coordinate. 

In Total Lagrangian Formulation Green strain tensor is used as a measure of strain and 
corresponding stress is Second Piola-Kirchoff stress (PK2) but we will use nominal stress 
also where required, the reason for this will be explained latter. A weak form of 
momentum equation is developed, which is known as principle of virtual work. 
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4.2.1 Some terms and designation: 

D is the deformation in terms of velocity or displacement; velocity is denoted by v(X,r), 
Cauchy stress is denoted by cr(X,r), rate of deformation is denoted by T)(X,t) and density 
is denoted by p(X,?),all these variable are dependent variables. Nodal coordinates in 
current configuration are denoted by , / = 1 to n (n is the number of nodes). Lower case 
subscript are used for component, upper case subscript are for nodal values. Nodal 
coordinates in deformed configuration are X ^, . 

In finite element method, motion x(X,/) is approximated by the relation 

x(X, 0 = N, (X)x,, (r) or x(X, t) = N, (X)X; (0 . (1) 

Where N, (X) is the interpolation (shape) function and Xj is the position vector of node I. 

We define the nodal displacement at the nodes by 

= or u,(t) = x,(t)-Xi . (2) 

The velocities are obtained by taking the material time derivative of the displacement, 
given by 

V, (X, t) = . ^4— - = (0^/ (^) or v{X,t) = u,{t)N,{X) . (3) 

at 

f'!"' and are the internal and external nodal forces. Define matrix B as 



dN, 

dXj 



4.2.2 Corotational Formulation: 


( 4 ) 


In structural elements such as bars, beams and shells fixed coordinate systems makes the 
analysis of structure difficult, for example a rotating rod. Initially cr^ is the only non-zero 

stress and vanishes, as the rod rotates its become awkward to express state of uniaxial 

stress in terms of global component of stress tensor. To overcome this difficulty 
coordinate system is embedded in the bar and is rotated with the rod, this is the concept 
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of corotational formulation and the coordinate system is known as corotational 
coordinate. This concept of corotational formulation is used in our F.E.A. 

4.2.3 Total Lagrangian Formulation: 

Governing Equations: 

In Total Lagrangian Formulation nominal stress P is used in momentum equation 
because the resulting momentum equation and its weak form are simpler than that for 
PK2 stress. PK2 stress is used in constitutive equation as P is not symmetric. 
Deformation tensor F does not vanish in rigid body rotation hence is not suitable for 
measuring strain in constitutive equation and constitutive equation is written in terms of 
Green strain tensor E. Consider a mass with domain Q and boundary defined byP, as 
shown in Figure 1 . 



Figure 14 Deformed and undeformed body 


Governing equations for Total Lagrangian formulation are as follows (subscript ‘0’ 
signifies variable with respect to undeformed position, while ‘inf, ‘exf and ‘kin’ 
superscripts denote internal, external and inertial values of the variables): 

Conservation of mass: 

pj = PqJq — Po • 

Here p is the density and J is the determinant of Jacobian between spatial and material 
coordinates. 
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Conservation of linear momentum: 


Vo-P + /?ob = Po«. (6) 

Here b is the body force. 

Conservation of angular momentum: 

F.P = P^-.F^ (7) 

Conservation of energy: 


PoM^"‘=F^P-Vo-q4-po^. 


( 8 ) 


w'"' is the rate of internal work done, q is heat flux, 5 is the power supplied by heat 
sources, here q = JF"' • q . 

Constitutive equation: 

S=S (E...etc.). P = S-F''. (9) 

Here S is PK2. 

Measure of strain: 


E = |(F'’-F-I). 

/ is identity matrix. 

Boundary conditions: 

n°Pj, = or e,. • n° • P = e- on F® , 


( 10 ) 


( 11 ) 


Here n is the unit normal, t is the surface traction, is the base vectors of coordinates. 


n, = m", our® , r; U r®„, = r®, rj n r° = o, for i=i to 


where n . is the number of space dimensions. 


Initial conditions: 


( 12 ) 
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P(X,0) = P„(X). 


(13) 


u(X,0) = uo(X). (14) 

Internal continuity condition: 

n]P,=QonTl. (15) 

4.2.4 Total Lagrangian Weak Form: 

We will develop weak form from strong form. Here the strong form consists of 
momentum equation (6), traction boundary conditions (12) and internal continuity 
condition (15). Defining spaces for test and trail functions as 

Sxi{X)eU„ u(W,/)eD (16) 


U is the space of kinematically admissible displacements Uq is the same space with the 
additional requirement that the displacements vanish on displacement boundaries. 

To develop weak form momentum equation is multiplied by test function and is 
integrated over the initial configuration and we get 








kfQo=0. 


Finally we get weak form as 


( Ju, u) - ( Ju, u) + (^u, u) = 0 V (Ju G D, 


(17) 


(18) 


Here we have 


= {SF^: FdQa • 

^.Vf/ « ~ 


(19) 

( 20 ) 


Sw’’’" = f Su-p^ndQ^ 


( 21 ) 
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From weak form (17) it is clear that work done due to internal force is equal to sum of 
work done by external force and inertial force, hence for finite element analysis of the 
model we can just obtain the value of internal force. 

Internal nodal force: 

From(19) ■J""' = • (22) 

From (4) and (22) 

• = (23) 

4.2.5 Linearization: 

Linearization of internal nodal forces is done and we derive the tangent stiffiiess matrix 
AT'"' for Lagrangian elements. Linearization of the constitutive equation is carried out 
with the continuum tangent moduli, the resulting material tangent stiffiiess matrix is 
called the tangent stiffiiess matrix. The tangent stiffiiess matrix is expressed in terms of 
C ^''- , the tangent moduli relating the rate of PK2 stress to rate of Green strain in the total 
Lagrangian framework. Tangent stiffness matrix is developed by relating the rates of the 

internal nodal forces to the nodal velocities J . Taking the material time derivative of 
(23) we get 



To obtain the stiffness matrix K'"^ strain rate P is expressed in nodal velocities using the 
constitutive equation and the strain measure, here we use material time derivative of the 
PK2 stress. Time derivative of the transformation F = S- F' is given by 

P=S-F"+S-F^ or Pi, =5-F;+S,,FJ . (25) 

fi-om (24) and (25) 
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(26) 

From (26) it is evident that the rate of internal nodal forces consists of two parts: 


1) The first part which involves the rate of stress (S) depends upon material 
response and leads to material stiffness matrix denoted by K"'°‘ . 

2) The second part involves the current state of stress S, and leads to geometric 
effects of deformation (including rotation and stretching). This term is called 
geometric stiffness and is denoted 

The internal nodal force can now be written as: 

. int mat . 

{ =f +f or f;r=f:;^‘+fr . in) 

First term in (27) material rate of internal nodal forces and second term is geometric rate 
of internal nodal forces. 

Material tangent stiffness: 

Voigt form of notation is convenient in developing tangent stiffness matrix as the tensor 
of tangent moduli is forth order tensor and cannot be handled readily by standard 

matrix operation. From equation (26) and (27) we get 

{ =J'Bj{S}dno (30) 

where {5} is rate of Pk2 stress in Voigt column matrix form. Now constitutive equation 
is written as 


Sy=C%E, or {5} = [C^"]{3E} (31) 

and the rate of the Green strain to nodal velocities in Voigt notation is given 
by {E} = Bo . From (30) and (3 1) we get 
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( 32 ) 


. int 

Finally the material stiffness matrix is obtained as 

K“ = |^B;[c“]B.dQ. or K^' = {^Bj, [C'-‘']B„rfn. . (33) 

Geometric stiffness: 

The geometric stiffness is obtained from relation = dN, /dX^ and from (26) and (27) 

Kr=lf (34) 

From (34) we see that PK2 stress is in the form of tensor, i.e. a square matrix. Each 
submatrix of the geometric stiffiiess matrix is unit matrix, hence the geometric stiffness 

A SCO 

matrix is invariant with rotation, i.e.K/j =KfJ", here the superposed denotes the 
geometric stiffness matrix on rotated coordinate system (corotational form). 

Material and geometric tangent stiffness matrix for an element: 

To obtain the material and geometric stiffness matrix for an element consider the 
element as a two-node rod element in two dimensions as shown in the Figure 2. The 
element is in uniaxial state of stress with the only non-zero stress along the axis of bar. 
To simplify the formulation, consider material coordinate system to coincide with the 
axis of the rod (Corotational formulation), with the origin of the material coordinate at 
node 1. The parent material coordinate is^ € [0,1] . The material coordinate is related to 

the element coordinates by 

X = (35) 

/ 

Here is the initial length of the element. The motion is given in terms of the element 
coordinates by 

fxl fi, x(f,f) = x,(/)N,(f). (36) 

bJ bi T2JK J 


where 
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( 37 ) 


1_X X 

K h 

y 

X 


1 

Initial Current 

Figure 15 Rod element in two dimensions in total Lagrangian formulation 
The Bq matrix from (4) is defined as 

[4,]-[aA',/3X,f=i[-l +1], (38) 

Let US take the coordinates of node 1 and node 2 as [x, ^ and [xj 3^2 ] respectively. 

Define Xj, = Xj - x, and j/ji = 3^2 " 9 = x^Jl an^smO = y^J I . 

Here ‘F is the deformed length of the element. In corotational formulation the only non- 
zero stress is along the axis of the rod. The X axis coincides with the axis of the rod, so 
S'], is always the stress component along the axis of the rod and in Voigt notation only 

the first row of the Bg matrix need to be considered. The By matrix is written as in Voigt 

notation as: 






34 


( 39 ) 



dN, dx 

'WW 

dN I dx 

dN, dx dN, dx 
dX dY'^'WM 


dX dX 
dN, dy 

dN, ^ dN, dy 
dX dY 


Now we havex,X - Xj, //() -/// q cos^, = >< 2 , //q = ///q sin0 . The first row of the 
Bq from (39) is given by 


B/=[A^/.x^>x ^,,xy^x] = [^,,xhi (40) 

From (37) we have 

Bo=[^B° Bj J = -^[-cos(9 -sin^ cos^' sin^]. (41) 

Material tangent stiffness matrix: 

From (33) and (41) we get 


K" 


AC 


si; ( I 


/. 


\hj 


cos^ 6 cos ^ sin ^ 

sin^ 9 


-cos^ 6 
-cos^sin^ 
cos^ 9 


symmetric 

Here is the undeformed cross section area of the rod element. 

Geometric tangent stiffness matrix: 

From (34) and (38) we get 


K«“'=- 



■+1 

0 

-1 

0 


0 

+1 

0 

-1 

io 

-1 

0 

+1 

0 


_ 0 

-1 

0 

+1_ 

holds 

for 

any 

orientation 


- cos 6* sin 6* 
-sin^ 9 
cos ^ sin ^ 
sin^ 9 


(42) 


(43) 


matrix is the sum of the material and geometric stiffness. 

In Updated Lagrangian form material and geometric stiffness matrix can be 


written as 
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cos^ 0 

cos 0 sin 0 

-cos^ 9 

-cos^sin^ 



sin^^ 

-cos^sin^ 

-sin^ 9 

/ 



cos^ 9 

cos ^ sin 0 


symmetric 



sin^^ 


+1 0 




I 


0 

-1 

0 


+1 

0 

-1 


-1 0 

0 -1 
+1 0 

0 +1 


( 44 ) 


(45) 


is the tangent moduli relating the Truesdell rate of the Cauchy stress to the rate-of 
deformation in the updated Lagrangian formulation, is the only Cauchy stress 
present. 


According to classical theory of rubber elasticity for non-Gaussian chain force on a chain 
strand having a vector length r is given by [4] 


1 . 




(46) 


Where N is the number of statistical chain segments /, is there length. is the inverse 
Langevin function. Now stiffness of this strand is obtained from df ! dr using the 
approximation for Langevin function 

Z’' ( 5 ) = 5(3 - 5^ ) / 1 - (5 is any variable) 

as 


WHl + Nl/ 


(47) 


Also there is Van der Waals force present in the model and the force for particular Van 
der Waals element is obtained as 


f.j = Ks/a 


(48) 
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where f ,j is the force in the Van der Waals element, K is the stiffness and s is the strain 
in the Van der Waals element. Following [4] 

a = (nf3)N''^. (Aq\ 


Final stiffness of an element is due to both strand stiffness and also Van der Waals 
stiffness. Hence (45) and (46) becomes 



9 

cos ^ sin ^ 

- cos^ 9 

-cos 9 sin 6* 

K, + K 


sin^ 9 

-cos (9 sin ^ 

-sin^ 9 

1 



cos^(9 

cos ^ sin 9 


symmetric 



sin^ 9 


(50) 


K 


am 


+1 0 


/ + /v 


vd 


0 

-1 

0 


+1 

0 

-1 


-1 

0 

+1 

0 


0 

-1 

0 

+1 


(51) 


We can now analyze the model of the polymer using (50) and (51) to get the strain 
energy vs. strain and stress vs. strain curves. 
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CHAPTER 5: RESULTS AND DISCUSSION 


5.1 Introduction and objective of the study; 

Polymer model developed in chapter 3 was used and the results obtained from these 
analyses were used to predict the mechanical behavior of amorphous polymers. The 
purpose of this work is to demonstrate that the model based on rubber elasticity has the 
potential of qualitatively and quantitatively predicting the properties of amorphous, 
glassy polymers. 

5.2 Actual stress-strain curve of an amorphous polymer: 

Figure 1 shows that there are three distinct regions on the true stress-strain response of a 
polymer. 

1) Initially the stress rises in an approximately linear manner as the applied strain 
increases. In this region the load bearing capacity of the polymer derives from 
both Van der Waals as well as molecular bonds in chain strands. 

2) At the yield point there is a abrupt fall in true stress. The fall in true stress is due 
to breakage to Van der Waals bonds which are weaker than strands and break 
earlier than strands. 



Figure 16 Actual Stress-Strain curve of a polymer showing different stages of deformation 
After breakage of Van der Waals bonds the entire load is taken by the strands. 
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3) Finally at large extensions the slope of the stress-strain curve increases again, i.e. 
a strain hardening effect occurs. This is because chains orient in the direction of stretch. 

5.3 Defining the problem: 

We had taken a regular grid of five hundred nodes (entanglement points) with periodic 
boundary conditions applied in x-direction. Polymer chains were obtained on this grid 
using S.A.W. Each strand is taken as a finite element in our model. Later on Van der 
Waals forces are applied between all nodes which add up to the stiffness of the model. 

5.3.1 Application of boundary conditions: 

Displacement boundary conditions are applied to the specimen in the y-direction. The 
nodes at the top of the specimen are given a displacement in the y-direction and the nodes 
at the bottom of the specimen are fixed and cannot move in y-direction. However, they 
are free to move in x-direction. The nodes at the extreme left of the grid are fixed in x- 
direction and the nodes at the extreme right of the grid are given a displacement in x- 
direction, so that the volume of the specimen remains constant. We had taken X as the 
stretch of specimen in the y-direction, as the two-dimensional model system is assumed 
to have a constant three-dimensional volume then if the specimen is stretched by X in y- 
direction then it is contracted by a factor of X~°’^ in the x-direction. Let n^be the 

displacement in y-direction and Zj the specimen length. If u is the rate of deformation 
per unit time in y-direction of the specimen then the incremental stretch dX^ in y- 

direction is given by u At , where At is the incremental time. Incremental stretch dX^ in 
x-direction is given by dX^/lX^^'^ , here X^ is the stretch in y-direction and is given by 
relation I + u^/Lq. This ensures that the deformation is volume preserving. 

5.4 Defining Strain energy: 

The area under the stress-strain curve up to a given value of strain is defined as the strain 
energy per unit volume consumed by the material in straining it to that value. This is 
easily shown by the relation. 
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( 1 ) 



where V is the volume of specimen, U is the strain energy, a and s are the uniaxial 
stress and strain respectively. P and L are the load and length of the specimen 
respectively. The tensile engineering stress in the simulated microscopic network, 
keeping the volume constant and keeping the assumptions of rubber elasticity is given by 
the relation. 


1 

cr = — 
V 




(2) 


The derivative should be carried out at energy minimum with respect to the coordinates 
of all the internal nodes at each incremental deformation step. 


5.5 Van der Waals bonds and their breakage: 

We had discussed about Van der Waals forces in the previous chapter. Here we will 
discuss the phenomenon of the breakage of these Van der Waals bond and show some 
results for different values of the stiffness of these Van der Waals bonds. The model is 
based on the Eyring activation rate theory [3], which takes into account the role of weak 

Van der Waals forces between chains. The network is elongated at a constant rate s and 
at temperature T along the y-axis. This causes straining of weak Van der Waals bond 
which break according to the kinetic theory of fracture [5], at the rate 

v = Texp[-{U-/3cT)/kT]. (3) 

Here t is the thermal vibration frequency and U and p are the activation energy and 

volume respectively, cr is the local stress and satisfies the relation a=—Yis , where s 

a 

is the local strain and K is the elastic constant of the bond, k is the Boltzman’s constant 
and T is the deformation temperature. The quantity a is given by the relation [1,2] 

a = (n/3)y. (4) 

Here N is the number of statistical chain segment and / is their length, n is the number of 
chain strands per unit volume. 
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As the Van der Waals bond break the external load is transferred to now effectively 
isolated chain strand. The simulation of the Van der Waals bond breaking is performed 
with the help of a Monte Carlo lottery that breaks a bond i according to the probability 

(5) 

^max is the rate of breakage of the most strained bond in the array. After each visit of 
a bond, the time t is incremented by 1/ [v^n(r)], where n(t) is the total number of intact 

bonds at time t. For analysis [1,2] r was taken to belO’^/ 5 . We have chosen a 
displacement increment of 10“*® m per time step. Deformation temperature (T) is taken to 

be 300 K. We have chosen [1,2] OlOO kJ/mol, andy0 = (26^)^ « = 3.04xl0“/m^and 
N= 5. Results for different values ofX are shown in the figures below. 



X (Angstrom) X (Aigstrom) 


(a) 


(b) 


Figure 1 Deformed polymer for only Van der Waals bond and K=lMPa 



XCAngslram) X (Angstrom) 

(a) (b) 

Figure 3 Deformed meshes for only Van der Waals bond and K=4 MPa 
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(a) (b) 

Figure 4 Deformed meshes for only Van der Waals bond and K=10 MPa 

Figure 2 shows deformed shapes of the polymer model with only Van der Waals bond, 
Stiffness of the bonds is taken equal to 1 MPa. The Figure 2(a) shows the shape of 
the polymer when the breakage of Van der Waals bond just starts (80% strain) and the 
Figure 2(b) shows the sh^e of polymer when some of the bonds are broken (90% strmn). 
The undeformed mesh is shown by dots in the background. 

Figure 3 shows the deformed meshes of polymer model with only Van der Waals bonds 
and stiffness (K) of bond equal to 4 MPa. Figure 3(a) shows the deformed mesh when the 
Van der Waals bonds just start breaking at a much lower strain than the case when K=1 
MPa (25% strain) and Figure 3(b) shows the deformed mesh when some of the bonds are 
broken (27% strain). 

Figure 4 shows deformed meshes for polymer model with only Van der Waals bond and 
stiffness (JQ of bonds equal to 10 MPa. Figure 4(a) shows the deformed mesh when 
bonds just start to break and Figure 4(b) shows the deformed mesh when some of the 
bonds are broken. For above value of stiffiiess the entire bonds break nearly at the same 
time as due to high stiffness the model behaves like a brittle material and fails without 
much deformation. As the stress-strain relation is linear for Van der Waals bonds the 
model with only Van der Waals bond deform in a regular manner. 
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Strain 

(C) 

Figure 5 Strain Energy-Strain curve for meshes with only Van der Waals bond and K=1 MPa, K=4 
MPa and K=10 MPa respectively 

Figure 5 shows the Strain Energy-Strain curve of polymer with only Van der Waals 
bonds and for stiffiiess values K=1 MPa [Figure 5(a)], K=4 MPa [Figure 5(b)], K=10 
MPa [Figure 5(c)]. From 5(a) it can be seen that after 95% strain the curve flattens out 
indicating that all the Van der Waals bonds are broken. This flattening of Strain Energy- 
Strain curve is obtained after 30% and 10% strain for K==4 MPa and K=10 MPa 
respectively. As the stiffiiess of the Van der Waals bonds is increased the model becomes 
stiffer and the bonds break at lower strain. This behavior is shown by the curves obtained 
for different values of stiffiiess of Van der Waals bond. 
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(c) 

Figure 6 Stress-Strain curve for meshes with only Van der Waals bond and K=1 MPa, K=4 MPa and 
K=10 MPa respectively 

Figure 6 shows the Stress-Strain curve of polymer with only Van der Waals bonds and 
for stiffness values K=1 MPa [Figure 6(a)], K=4 MPa [Figure 6(b)], K— 10 MPa [Figure 
6(c)]. The abrupt drop of stress to zero at particular value of strain for various stress- 
strain curves indicates that at this point entire Van der Waals bonds are broken. As the 
value of stififiiess increases the slope of stress-strain curve also increases. It is clear from 
the jcurves that the Van der Waals bonds with less stiffness value can sustain larger 
deformation than the bonds that are more stiff. 
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5 .6 Strands and their breakage; 

As discussed in chapter 3 and chapter 4 strands that are stretched uniaxialy follow the 
relation 

kT 

f = -jL-\rlNl). (6) 

Here/ is the force in a strand; / is the length of statistical chain segment is taken to be 

10 A [1,2]. For variable x we have approximated the Inverse Langevin function by the 
relation!, '(x) = x(3 -x^) /(I- x^). Clearly as x->l the function !“'(x) ->Qo. Plots for 
N= 5, 12, 20 are shown below. Here no Van der Waals bonds are considered, mesh 
consists only of strands. Since the stress function depends upon inverse Langevin 
Function hence the strain energy-strain curve shows a steep increase with strain. The 
stress-strain relation is also non-linear. Polymer model is deformed to 30% and 60% 
strain and are the deformed meshes are plotted. The final strain energy-strain curve and 
stress-strain curve is obtained for 60% strain. At this value of strain the most of the 
molecular bonds are broken; we break a strand when the draw ratio of the strand 

approaches [19] as at this value of draw ratio the stress approaches infinity 
(cr->oo). 

Figures 7(a & b), 8(a & b) and 9(a & b) show the deformed meshes for polymer model 
with only strands and number of statistical chain segments (N) equal to 5,12 and 20 
respectively. The two deformed meshes for each figure are plotted for 30% and 60% 
strain respectively. At 30% strain some of the strands are broken and later on most of the 
strands break giving rise to material failure. The strand is assumed to be broken when its 

draw ratio approaches a value ofVv . As can be seen fi'om deformed meshes that as 
value of N increases less number of strands are broken for same strain value this is due to 
the fact that at these strain value the strands having higher value of N don’t attain the 

draw ratio of y[N . 
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(C) 

Figure 10 Strain Energy-Strain curve for meshes with only strands and N=5, N=12, N=20 
respectively 


Figure 10(a, b, c) shows the Strain Energy-Strain curve for polymer model with only 
strands and number of statistical chain segments (N) equal to 5,12 and 20 respectively. 
Strain Energy is a function of Inverse Langevin function hence the curve of strain energy- 
strain is expected to be as obtained and hence the curve is non-linear in nature. The 
polymer models were deformed to 60% strain to get these plots. 
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Strain 


(c) 

Figure 11 Stress-Strain curve for meshes with only strands and N=5, N=12, N=20 respectively 

Figure 10(a, b, c) shows the Strain Energy-Strain curve for polymer model with only 
strands and number of statistical chain segments {N) equal to 5,12 and 20 respectively. 
The Stress-Strain curve is non-linear in nature and as the number of statistical chain 
segments are increased from 5 to 20 the polymer model attains less maximum stress for 
same strain value also the non-linearity of the stress-strain curve decreases with the 
increase inA^. 
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5.7 Actual polymer model with both Strands and Van der Waals forces: 

Now results are obtained for both strands and Van der Waals bonds present in the model. 
Due to Van der Waals force we get the first portion (1) of stress-strain curve in Figure 1 
as nearly linear then as the Van der Waals bonds start breaking we get the softening part 
(2) of the curve and finally for large strain we get the strain hardening part (3). In these 
following analysis K=4 MPa is kept constant and N is varied. 



Strain Strain 


Figure 12 Strain energy-strain and stress-strain curve for N=12 and K=4MPa 

Figure 12 shows the strain energy-strain curve and stress-strain curves for polymer model 
with both Van der Waals bonds and molecular strands present with N=12 and stiffness 
(K) equal to 4 MPa. From the stress-strain curve it can be seen that up to aroimd 20% 
strain Van der Waals bonds do not break and at 27% strain all Van der Waals bonds are 
broken. After breaking of Van der Waals bonds strain hardening of the polymer takes 
place. 


X lo ’* 



Figure 13 Strain energy-strain and stress-strain curve for N=20 and K-4MPa 
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Figure 13 shows strain energy-strain and stress-strain curves for polymers model with 
both Van der Waals bonds and molecular strands present with N=20 and stifl&iess (K) 
equal to 20 MPa. Van der Waals bonds start breaking at around 25% strain and all the 
bonds are broken at 30% strain. Stress at same value of strain for V=20 is less than stress 
for V=12 which is expected behavior following the curves which were obtained only for 
strands. 



X (Angstrom) X (Angstrom) 


Figure 14 Deformed meshes when all the Van der Waals bonds are present and when all the Van der 
Waals bonds are broken for N=5 and K==4 MPa 

Figure 14 shows the deformed meshes for polymer model with Van der Waals bonds 
and strands present. N=5 and K=4 MPa is taken, the left hand side figure shows the 
deformed mesh when all the Van der Waals bonds are present (25% strain) and the right 
hand side figure shows the deformed mesh when all the Van der Waals bonds are broken 
(30% strain). It can be seen in Figure 16 that due to presence of Van der Waals bonds the 
polymer model deforms in a regular fashion and as the Van der Waals bonds are broken 
the deformation of polymer is in an irregular fashion as can be seen fi’om Figure 17. In 
Figure 15 the polymer is deformed up to 60% strain. In this figure it can be seen that 
most of the strands are broken. 
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Figure 15 Deformed mesh after 60% strain for N=5 and K=4 MPa 




Figure 16 Strain energy-strain and stress-strain curve for N=5 and K— 4MPa 

Figure 16 shows the strain energy-strain and stress-strain curves for polymer with Van 
der Waals bonds and strands present. N=5 and .^=4MPa is taken for the analysis of the 
polymer model. 


fWTfgr 2)0 A- 
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Figure 17 Plot showing Stress-Strain, Strain Energy-Strain, Number of Van der Waals bonds 
broken-Strain for N=5 and K=4 MPa 

Figure 17 shows the plot for Stress-Strain, Strain Energy-Strain, Number of Van der 
Waals bonds broken-Strain for N=5 and K==4 MPa. From the plot it is clear that the strain 
hardening starts when all the van der Waals bonds are broken here also the curve of strain 
energy-strain curve slope changes. 
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CHAPTER6: CONCLUSION 


A model is proposed for deformation of glassy amorphous polymers and the analyze of 
the same is done using Finite Element Analysis. The parametric data is taken from the 
reference papers and the results qualitatively mimic the actual behavior of amorphous 
polymers. 

6. 1 Scope for future work 

As the data for analysis of the polymer model (for example the data taken in the equation 
of the breakage of Van der bonds using kinetic theory of fracture) is taken from reference 
papers attempt should be made to verify data from experimentation. Moreover atomistic 
analyze can be used for better understanding of the Van der Waals bond breakage 
phenomenon. 

We have done a qualitative analysis of the polymer model which can predict the nature 
of the stress-strain behavior of the actual model, but the values of stresses at a particular 
strain value do not match with the actual stress value of a polymer under study. Hence the 
model should be extended to do a quantitative analysis. 
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